Are High-Redshift Quasars Blurry? 
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ABSTRACT 



It has been suggested that the fuzzy nature of spacetime at the Planck scale 



may cause lightwaves to lose phase coherence, and if severe enough this could 

blur images of distant point-like sources sufficiently that they do not form an 

Airy pattern at the focal plane of a telescope. Blurring this dramatic has already 

been observationally ruled out by images from Hubble Space Telescope (HST), but 

I show that the underlying phenomenon could still be stronger than previously 

considered. It is harder to detect, which may explain why it has gone unseen. A 

systematic search is made in archival HST images of among the highest known 

redshift quasars. Planck-scale induced blurring may be evident, but this could 

■^j- , be confused with partially resolved sources. 

O 

\Q ', Subject headings: time — gravitation — quasars: general 

O 



1. Introduction 

Characterizing the microscopic properties of spacetime using images of distant sources 
was first proposed in a Letter by Lieu & Hillman (2003). This compelling possibility is a 
consequence of a phenomenon likely to be found at the Planck scale, where lengths shrink 
to Zp ~ 10~ 35 m and time intervals diminish to tp ~ 10~ 44 s. Following Ng, Christiansen, 
& van Dam (2003), here distance measurements / are uncertain by 51 > Zp(Z//p) 1_a (with 
similar uncertainties for time), where the parameter a specifies different quantum gravity 
models. A natural choice for a is 1, but ^ is also possible (Amelino-Camelina 2000), and | is 
preferred (Ng 2002), because it is consistent with the holographic principle of t'Hooft (1993) 
and established black-hole theory (Beckenstein 1973; Hawking 1975). An electromagnetic 
wave travelling a distance L from source to observer would be continually subjected to these 
random spacetime fluctuations, which following Ng, Christiansen, & van Dam (2003) leads 
to a cumulative statistical phase dispersion of 

A0 O = 2 7 ra |L 1 - a , (1) 
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where A is the observed wavelength, and ao is close to unity. This presents a means of 
detection in an image of a point source. Ragazzoni, Turatto, & Gaessler (2003) point out 
that A0 O can be interpreted as an apparent angular broadening of the source, and over a 
sufficient distance might grow enough to obliterate the diffraction pattern at the focal plane 
of a telescope. 

Because this has not been seen in Hubble Space Telescope (HST) data, Lieu et al. and 
other authors (Ragazzoni, Turatto, & Gaessler 2003; Ng, Christiansen, & van Dam 2003) 
have concluded that the effect, if present, is too weak to be observed. In Section 2 I expand 
on previous work to account for the shorter wavelength of emitted light, which predicts 
stronger blurring. But observing this is still problematic, because it must be disentangled 
from a partially resolved source. A test is outlined in Section 3 which makes use of the best 
available archival data, a HST snapshot survey of z > 4 quasars. These data are discussed 
in Section 4, after which a new analysis in the context of blurring by Planck-scale effects is 
presented in Section 5. 



2. Limits on Blurring Induced by the Planck Scale 



Equation 1 assumes that the emission and detection wavelengths are the same, and 
so represents a minimal estimate for blurring. For cosmological distances, photons left the 
source at the shorter wavelength for which the blurring action of the Planck scale should 
be stronger. Consider the following: 

A0 max = 2na (l + z)%L l ~ a = (1 + ^)A0 O - 



(2) 



One might call this a maximal estimate, because it assumes the bluest possible photons 
propagate a distance L and are then detected at wavelength A. It is instructive to write 
equation 2 in a different form. By taking its differential with respect to z, it is possible to 
rewrite it as 
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= A0 Lo . s . + A<j) z . 
where the luminosity distance L is given by 



q z - (1 - q )(\/l + 2q z - 1) 



(3) 
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is the deceleration parameter. By comparison with equation 1, the 



first integral can be recognized as the observed phase dispersion of photons of wavelength A 



- 3- 



arriving from all redshifts up to the source. In practice these could conceivably come from 
intervening sources along the line of sight, hence the label A(f>\. 0mBm . The second integral must 
then be the remaining phase dispersion A<p z associated exclusively with photons redshifted 
to the observer's wavelength. The latter scenario seems the more plausible. This is plotted 
in Figure 1 for az = 4 point source as a function of observed wavelength using the standard 
LCDM cosmology (Qa = 0.7, = 0.3, and H = 70 km s _1 Mpc -1 ) assumed throughout 
this paper. Three other possibilities are also plotted: minimal blurring A0 O , the combination 
of minimal blurring of the source plus that of intermediate sources along its line of sight 
A0 O + A0i. o . s ., and maximal blurring A0 max , for a total of four. Following Ng, Christiansen, 
& van Dam (2003) I use a = | and a — 1. Reducing (increasing) a shifts these lines in 
parrallel towards lower left (upper right) in this plot, as does increasing (decreasing) a. The 
most dramatic changes come with variation in a, by virtue of it appearing as an exponent 
in equations 1 and 3. 

Can blurring be constrained by current observations? Figure 1 illustrates at what wave- 
lengths an answer may be most convincing. The highest spatial resolutions of some current 
telescopes are indicated. Chandra routinely observes sources to a resolution limit better than 
that imposed by either equations 1 or 3 for a — |. But this conflict is easily avoided, by 
either suitably increasing a (0.7 would be sufficient) or reducing a - or both. This flexibility 
is allowed because current telescopes operating shortward of the optical (including Chandra) 
do not form diffraction-limited images, and so it may be more difficult to distinguish blur- 
ring induced by the Planck scale from simply an unresolved source such as a galaxy. The 
dashed lines indicate the angular size A6 of extended z = 4 sources for a range of physical 
diameters. Thus, it may be best to look for blurring in the optical. HST probes a regime 
predicted to be blurry, which is indicated by the shaded regions. It can also resolve sources 
confirmed to be physically small at other wavelengths. For example, if a = | and subject to 
the maximum blurring limit, no z = 4 source should appear sharper than about C'l across 
at 800 nm, broad enough to affect the HST diffraction pattern. An object with z = 4 would 
need to be 500 pc across to mimic this level of broadening. This is smaller than a typical 
galaxy, but not quasars. Even if it were a factor of about 2.5 less, HST could still marginally 
discriminate between Planck-scale-induced blurring and the effects of a 200 pc source. But 
a ~ 25-m space telescope would be needed to rule out the minimum set by equation 1, by 
resolving objects as small as 50 pc. 



3. Observational Plan and Sample Selection 

Ragazzoni, Turatto, & Gaessler (2003) have already confirmed that the most severe 
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blurring discussed in Section 2 is not seen in HST images of a z = 5.4 Hubble Deep Field 
galaxy. It is worth considering how weak blurring could be and still be ruled out for HST. 
Ragazonni et al. point out that small perturbations of the phase should only cause a drop 
in image Strehl ratio S. This is the ratio of the image peak to that of the diffraction spike of 
an unabberated telescope. In this case, and where A0 is comparable to the telescope optical 
aberrations, the Marechal approximation to the Strehl ratio applies, and 



cxp 



Figure 2 is a plot of point-source Strehl ratio as a function of redshift for an unaberrated 
telescope with D = 2.4 m observing at A = 800 nm. The shaded regions have the same 
meaning as in Figure 1, again plotted for a — | and = 1. If the object is not a true point 
source but is actually extended, equations 1 and 3 do not strictly apply. In this case the 
Strehl ratio of the image would also drop due to the partial resolution of the source. For an 
object of angular size A6 close to the diffraction limit of the telescope one would expect a 
Strehl ratio of 

S « exp [ - (A0-J J (6) 

in the absense of Planck-scale effects. The dashed lines in Figure 2 represent equation 6 for 
various source diameters. 

Judging from Figure 2, to test if high-redshift sources are indeed blurred, one could 
look at a large homogeneous sample of compact high-redshift objects imaged with HST. In 
principle, the test is simple: Is there a downward trend in Strehl ratio with increasing source 
redshift? In practice, two factors must be well controlled. First, the sample must be of 
uniformly small sources of very high redshift. Figure 2 indicates that for any Planck-scale 
effect to be detected with HST the sample must extend to redshifts well beyond z — 4. 
Second, the point-spread function (PSF) must be well characterized because the Planck- 
scale-blurring signal is relative to the expectation for diffraction-limited imaging. 

As it turns out, there is already at least one archival HST dataset that meets these strict 
critera: the Advanced Camera for Surveys (ACS) High-Resolution Channel (HRC) snapshot 
survey of SDSS quasars. This was a program of imaging to detect lensed companions among 
a sample of 4 < z < 6.3 quasars. None were found despite the excellent spatial resolution 
afforded by ACS HRC. Because depth was not as important as spatial resolution, these 
images are relatively shallow, which should minimize confusion by any extended host galaxy. 
Also, Kaspi et al. (2005) suggest that none of the SDSS quasars should have a broad-line 
region larger than a few parsecs, which would make equation 6 negligibly less than unity 
for HST. Thus, this sample should satisfy the first condition. That they also satisfy the 
second is a property of both HST PSF stability and the Nyquist-sampled ACS HRC pixel 
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scale. Neither the Space Telescope Imaging Spectrograph (STIS) nor the Wide Field and 
Planetary Camera 2 (WFPC2) can fulfill this second constraint. 

4. Data and Reductions 

The archival HST ACS snapshot survey of Sloan Digital Sky Survey (SDSS) sources 
(Proposal: 9472, PI: Strauss) includes all SDSS quasars with z > 4; some of which have 
z > 6. They were obtained in snapshot mode through either the F775W filter (central 
wavelength 761 nm, 95 sources, 3.9 < z < 5.4) or F850LP (869 nm, 4 sources, 5.8 < z < 6.3). 
All were imaged with the HRC, which has a pixel scale of (X'0246. 

These data have already been discussed in detail in Richards et al. (2004, 2006), and 
those reductions are closely followed here. Pipeline-processed images were downloaded 
from the HST archive. This provides standard corrections for overscan, bias, and dark- 
current, flatfield division, bad-pixel masking, cosmic-ray removal, and photometric calibra- 
tion. The position of the quasar was determined to the nearest pixel using the IRAF 1 
program IMEXAM. Next, the Tiny Tim (v6.3) program (Krist 1995) was used to generate 
an appropriate PSF for each. Due to the strong positional dependence of the HRC field 
distortion this produces a more accurate PSF than can be determined by images of stars. 
Color dependence is also properly accounted for, by inputing a spectral energy distribution 
(SED) and convolving this with the HST filter curve. Synthetic PSFs were generated for 
each quasar using a redshifted template SED based on a composite SDSS QSO spectrum 
(Vanden Berk et al. 2001). By shifting in the Fourier domain, the sub-pixel position and 
peak of the PSF were allowed to float relative to the image, and minimized based on the 
sum of the square of residuals. The result is an unique PSF for each quasar, correct for its 
detector position, filter bandpass, and SED. 

All of the images were then co-aligned to the center of the nearest pixel. Figure 3 shows 
the result for SDSS J0836+0054 [z = 5.82), which will serve as an example. Each box is 3"x 
3" with the same orientation as in Richards et al. (2004) (their Figure 2). All images have 
the same grey-scale stretch. The non-circular first Airy ring is clearly evident, and further 
rings would be visible in a harder stretch. The diffraction pattern of a circularly symmetric 
2.4-m diameter pupil with a 0.7-m central obstruction is also shown. Note that this excludes 
any aberrations associated with the ACS camera, which explains its circular symmetry. The 
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same SED as the quasar was used. It is relative to this ideal diffraction pattern that - by 
definition - the Strehl ratio of the quasar and PSF is to be determined. Next to this is 
the PSF residual (the model residual will be discussed later, in Section 5). This residual 
is robust against variation in the template SED, both in relative line strengths and choice 
of continuum power law. It is also, as can be expected, comparable to that obtained by 
Richards et al. and reveals no obvious host galaxy or lensed component. But underlying 
effects due to microlensing cannot be ruled out, although it is not clear what effect they may 
have. Results are similar for the other sources. 

Strehl ratio was then measured, which is not reported by Richards et al. (2004, 2006). 
Each image and PSF was normalized to a flux of unity based on synthetic aperture pho- 
tometry, and the height of its peak determined relative to the ideal diffraction pattern. The 
accuracy of the final measurement is to within a few percent, which is the level of photometric 
uncertainty. 

5. Analysis and Results 

Figure 4 is a plot of Strehl ratio; crosses indicate the PSF. The unfortunate division 
between the F775W and F850LP data at z — 5.5 is evident. Even if one ignores this, a 
clear downward trend with increasing source redshift can be seen, despite some scatter in 
the data. Overplotted is a linear least-squares fit to these, which has slope of —0.05 ± O.Olz 
(l-o" errors). That some of the lowest redshift quasars have Strehl ratios higher than the 
PSF indicates that the limits of the telescope have been reached here. The differences of a 
few percent are consistent with the uncertainties in the measurements. Thus, the decrement 
below the PSF Strehl ratio for the z > 5 redshift sources is probably real. The situation for 
the z > 5.5 sources is more secure. Limits predicted by equations 1 and 3 assuming a — | 
and a = 1 are overplotted. The break is due to the change in filter effective wavelength, 
from 761 nm to 869 nm. For this choice of a and a, both A0 max and A0 O + A0i. o s . are 
clearly ruled out. But blurring associated with just redshifted sources is not. It is striking 
how closely the maximum Strehl ratio follows this curve. And of those quasars for which 
the PSF should be sharp enough to allow it, none are found in the region that it borders. 
If this is correct, it places a tight constraint on a. Even if ao = 10, a need only grow by 
2% (to 0.68) and still be in agreement with the data. Although also possible, this result is 
not easily explained by the size of the sources alone. It would seem that the intrisic sizes of 
these quasars (either their narrow-line or broad-line regions - or both) would need to occupy 
a very narrow range for this to happen, between about 200 pc and 300 pc. 

In a further simple test of the internal consistency of the z > 5.5 results, a simple 
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"blurred" quasar model was generated. Guided by the good fit of A0 2 in Figure 4, the 
existing PSF was cropped slightly below its peak, at a cutoff given by cutoff = peak x 

exp — ^A0 2 yj . A value A = 869 nm was used, for D = 2.4 m. This naive model is 
indeed a better fit to the observed Strehl ratio, indicated by the open triangles in Figure 4. 
The effects on the residual are benign, which can be seen in Figure 3. This is also evident in 
slices along the x and y axes, which are plotted in Figure 5. This confirms that the blurred 
quasar model is a better fit than the original PSF, accomplished without adversely affecting 
the residual. This is reassuring, as it demonstrates that these results are not in conflict with 
those of Richards et al. (2004). 

In summary, I have searched for blurring induced by the effects of the Planck scale in 
HST images of high-z quasars. Although blurring may be seen, if real, it is just at the 
observable threshold. The test here is more sensitive than previous ones, which only looked 
for compact sources lacking diffraction rings. The shorter wavelength of emitted light has 
also been accounted for, which if correct, restricts a to be 0.68 even if a is as big as 10. 

The implications of blurring are significant, so it is worth looking for. A true detection 
could point to a successful quantum gravity theory. But detection is elusive because the 
signal is confused with that of a partially resolved source. It is hoped that this work will 
encourage further tests for the effect with HST. A first step would be to re-observe the four 
z > 5.5 quasars in this sample with the HRC and the F775W filter. Observations of higher- z 
quasars (as they become known) could ultimately confirm or refute the current result. 
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Fig. 1. — Accumulated phase dispersions in arcseconds as a function of wavelength for a 
point source with z — 4, according to equations 1 and 3 with a — | and a — 1. By com- 
parison, the resolution limts of some telescopes are indicated: Chandra X-ray Observatory 
(Chandra), Roentgen Satellite (ROSAT), Far Ultraviolet Spectroscopic Explorer (FUSE), 
adaptive optics on an 8-m-class telescope (AO), Smithsonian Submillimeter Array (SMA), 
Very Large Array (VLA), and Very Long Baseline Interferometry (VLBI). The sizes of ex- 
tended objects of various physical diameters are shown as dashed lines. Shaded regions 
bordered by equations 1 and 3 (and two intermediate subregions formed by either excluding 
or combining their terms) are potentially precluded from observation. 
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Fig. 2. — Precluded Strehl ratio for blurred point sources (shaded regions) and predicted 
Strehl ratio for partially-resolved extended sources (dashed lines) as a function of redshift 
for HST observing at 800 nm. Shaded regions have the same meanings as in Figure 1. 
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Fig. 3. — Images of the quasar, PSF, ideal diffraction pattern, and residuals for one of the 
QSOs in the sample. Results are similar for the others. 
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Fig. 4. — A plot of Strehl ratio (open squares) for the sample. The crosses are the PSF 
for each quasar, giving an optically corrected Strehl ratio via the Tiny Tim software. There 
is a noticable trend towards lower observed Strehl ratio with increasing redshift. A linear 
least-squares fit to all data is overplotted. No correction for the change in filter at z — 5.5 
has been made. Models of Planck-scale induced blurring are shown; the shaded regions have 
the same meanings as in Figures 1 and 2. The break is due to the shift in filter central 
wavelength. Open triangles indicate the "blurred" quasar model of Section 5. The region 
bordered by A(p z is avoided, and so cannot be ruled out by the observations. 




Fig. 5. — Slices along the x and y axes of the same quasar shown in Figure 3 (dotted line). 
Overplotted are the PSF residual (thin solid line) and model residual (thick solid line). All 
are normalized to the peak intensity in the quasar image. The scale has been restricted to 
0.6 and below to better display the residuals. The model, although still not perfect, is a 
better fit to the quasar. 



